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1 Introduction 



Computing n-point functions of local composite gauge-invariant operators in = 4 supersymmetric 
Yang-Mills theory, SYM4, is an important problem because these n-point functions are among the 
basic objects on which the AdS/CFT correspondence can be tested ^ It is also a hard problem, even 
at tree-level, if only because of the combinatorial complexity of the operators involved. However, 
developments over the past few years, starting with [3] and subsequent works, raise the hope that the 
methods of classical and quantum integrability can be used to solve this problem, at least in the planar 
limit. For a comprehensive review of integrability in SYM4 and AdS/CFT, see [4] and references 
therein. For shorter review, see [5]. 

In this work, we focus on operators {O} that are composed of fundamental fields in the scalar 
sector of SYM4. Representing these fundamental fields as matrices in the adjoint representation of 
su{Nc), {O} are traces of products of Nc x Nc matrices. Further, in the planar limit that we are 
interested in, Nc — > 00, multi-trace operators are suppressed by factors of 1/Nc, and one can take 
{O} to be single-trace operators. 

The weak-coupling limit 

In weakly-coupled, perturbative Yang-Mills theory, the computation of the n-point functions is a well- 
defined problem. Following Okuyama and Tseng [6], it is sufficient at tree-level to count all possible 
planar sets of Wick contractions between the operators involved. Apart from normalization factors, 
the essential object in a 3-point function is the structure constant which, up to a normalisation is a 
tri-Unear form in the Hilbert space of states, which we call the cubic vertex, in analogy with string 
field theory. In [7], Roiban and Volovich showed that these n-point functions reduce to scalar products 
of spin-chain states constructed using the algebraic Bethe Ansatz [8]. A systematic study in the case 
of three operators that belong to (different) su{2) sectors was presented by Escobedo, Gromov, Sever 
and Vieira [9]. The tree-level correlation function of three su{2) operators was expressed in [9] in 
terms of scalar products of off-shell Bethe states^ of XXX spin- 5 chains. This method is known as 
"tailoring". Furthermore, it was shown in [10] that the 3-point function can be recast in terms of scalar 
products of an off-shell state and an on-shell state and thereby can be evaluated in determinant form. 
We refer to this method as "freezing". 

The semi-classical limit 

We are interested in computing the correlation functions of long operators that are dual to semi- 
classical string states in AdSs x S^. The semi-classical (heavy) operators {Ogc} are associated with 
classical solutions of the string cr-model [1 1-13]. For a review see [14]. In spin-chain terms, the oper- 
ators {Oi} are eigenstates of the spin-chain Hamiltonian. In other words, they are functions of rapidity 
variables that satisfy Bethe equations. For such operators, the Bethe roots condense into several cuts 
(macroscopic Bethe strings) in the complex rapidity plane [11]. The n-point functions of semiclas- 
sical operators {Osc} are particularly interesting, as they can be compared with the corresponding 
correlation functions computed on the string theory side. Computing n-point functions of {Osc} in 
the string theory was addressed in [15-22]. However, the only case when the complete answer is 

' One of the early tests of the conjecture was to check that the tree-level n-point functions of the BPS operators coincide 
with those in supergravity [1,2], 

^If the magnon rapidities satisfy the Bethe equations, the Bethe state is called on-shell, otherwise the Bethe state is called 
off-shell. 
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known is that of two heavy and one light operators [16,23,24]. The same configuration (heavy-heavy- 
light) was considered on the gauge theory side by Escobedo et al. [25,26], and in the sl{2) sector by 
Georgiou [27]. They used a coherent state approximation for the two heavy operators in the sm(3) 
sector. Comparison with the Frolov-Tseytlin limit [28] of the string theory result [16,23,24] showed 
a perfect match. 

The general case, when all three asymptotically-long operators are non-BPS, the complete answer 
for the three-point function is known only for weak coupling, and for special choice of the operators. 
In [29], Gromov, Sever and Vieira presented a thorough analysis of the case of one BPS and two 
non-BPS heavy fields from the 5-^(2) sector. In spin-chain terms, BPS operators are characterized 
by trivial Bethe roots that are set to infinity. The main result of [29] is an analytic contour integral 
derived from Korepin's sum expression for the scalar product of two off-shell states [30]. In [31, 
32], the determinant expression obtained in [10] was used to solve the problem in the general case 
of three heavy non-BPS operators. In [33-35], it was argued that this solution gives the two-loop 
approximation for the 3-point function as well. 

Outline of contents 

In Section 2, we classify the 3-point functions such that at least one operator is from the sn(3) sector. 
In 3, we recall the formulation of the su{2) 3-point function, including the cubic vertex, in determinant 
form. In 4, we use a result of [36] to generalize the determinant expression of the 3-point function, and 
the cubic vertex, obtained in [10], where all operators belong to scalar su{2) sub-sectors, to the case 
where two of the operators belong to su(3) sub-sectors while the third belongs to an su{2) sector. To 
use the result of [36], we need to take one set of Bethe roots of one of the su{^) operators to be trivial 
(sent to infinity). In 5, we recall the how the s'u(3) 3-point function was written in [3 1, 32, 37] in terms 
of certain functionals in order to be able to compute its semi-classical limit. In 6, we write the su(3) 
s-point function in terms of the same functionals as in the previous section. In 7, we compute the semi- 
classical limit of the su{'i) 3-point function of three non-BPS operators, under an assumption that 
allows us to compute the semi-classical limit of the norm of an sn(3) Bethe eigenstate. Appendix A 
contains a brief introduction to the nested coordinate Bethe Ansatz which is needed for the 'tailoring' 
approach to the su(3) 3-point function. Appendix B includes details of the 'tailoring' approach to the 
sm(3) 3-point functions. Appendix C discusses the properties of the functional forms that are needed 
to obtain the semi-classical limits. 

2 3-point functions with at least one su(3) operator 
2.1 The structure constant in A/* = 4 SYIVI 

The 2-point and the 3-point functions are determined, up to multiphcative factors, by conformal in- 

variance, 

mx.)0,{x,)) = (2.1) 
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where Nc is the number of colors, A is the 't Hooft coupling and the three factors Mi depend on 
the normalization of the operators {Oi{xi)}. The structure constant Ci23(A) does not depend on the 
normalization. 

The factors Li, equal to the number of fundamental operators in Oi, account for the cyclic rotations 
of the trace operator Oi. The structure constant Ci23(A) has the perturbative expansion 

Ci23(A) = cig + AcS + ... . (2.3) 

To compute the tree-level structure constant c[^^ using the method of [9], one needs the 1-loop 
wave functions. At 1-loop level, the operator Oi is represented by a Bethe eigenstate with energy Aj. 




Figure 1: Planar contractions contributing to the tree-level 3-point function. The con- 
tractions {XX), (YY) and (ZZ) are represented respectively by black solid lines, blue 
solid lines and dashed lines. 

In a su{2) sector, there is only one non-trivial configuration of 3-point functions. In the presence of 
one or more operators from an su{3) sector, the structure of the 3-point functions becomes richer and 
we need to classify the set of possible non-trivial configurations of structure constants. An example 
of a set of planar contractions is given in Fig. 1. The contractions ( XX ), ( YY ) and ( ZZ ) are 
represented respectively by solid lines, (red) wavy lines and dashed lines. 

Let us introduce some conventions. There are several possible choices of an su{3) sector, which 
correspond to a choice of three distinct complex scalar fields X, Y, Z, X, Y, Z, with pairs of mutually 
conjugated fields, like Z and Z, excluded. When only two types of non-conjugate scalar fields are 
chosen, the composite operator belongs to an su{2) sector. If Oi, O2 and O3 belong to su{a), 
su{(3) and 5^(7) sectors respectively, then the corresponding 3-point function of type {a, /3, 7}. By 
permutation invariance, the order of a, /3, 7 is irrelevant. 

We represent the different classes of correlation functions schematically by specifying the differ- 
ent types of Wick contractions between pairs of operators. For example, the correlation functions 
corresponding to Fig. 1 belong to the class {2, 3, 3} in Fig. 2. We call the operator at the bottom Oi, 
the one at right Os and the one at left O2. Exchanging a scalar field and its complex conjugate in all 
the operators does not change the value of the structure constant. This enables us to choose Oi such 
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Figure 2: Schematic representation 
of the type-{2, 3,3} correlation func- 
tion from Fig. 1. 



Figure 3: The two type-{3, 3, 3} structure con- 
stants. 




Figure 4: The remaining two type-{2, 3, 3} structure Figure 5: The type-{2, 2, 3} structure 

constants. constant. 

that it contains only the scalar fields X, Y and Z. Since we are interested in the large Nc limit, only 
planar contractions are retained. 

We start by classifying the type-{3, 3, 3} structure constants. In this case, there are two non- 
trivial inequivalent configurations, as is shown in Fig. 3. Deleting one line, that is, one type of Wick 
contractions, from each of these two configurations, one obtains type-{2,3,3} structure constants. 
There are three such configurations, as shown in Fig. 2 and Fig. 4. 

Deleting one line from the configurations in Fig. 4, one obtains a type-{2, 2, 3} or a type-{2, 2, 2} 
3-point functions. The latter is a pure su{2) 3-point function of the type studied in [9, 10, 31, 32]. 
There is one configuration of type-{2, 2, 3}, as in Fig. 5. To summarize, there are six non-trivial types 
of 3-point functions with at least one su{3) operator. 

2.2 Tailoring the tree-level structure constants 

Following [9], we construct the structure constant in three steps. 1. We split the algebraic Bethe 
Ansatz representation of each spin chain into two: a left sub-chain and a right sub-chain. 2. Consider- 
ing each spin chain to be an in-state, we "flip" each left sub-chain from an in-state to an out-state. 3. 
We take the scalar products of the left sub-chain state of Oi with the right sub-chain state of Oj+i mod 3 
(i = 1, 2, 3). Finally, we normalize the three external states. Further details on the tailoring procedure 
are in appendix B. We give our setup data in the table below. 



Operators 


Length 


Rapidities 


No. of Rapidities 


Partitions of Rapidities 


Oi 


Li 


Ul,U2 


#ui=A^i,#U2=Mi 


u'^ U = Ui, U2 U = U2 


O2 


L2 


Vl, V2 


# Vi=Ar2, # V2=M2 


v'^ U V{ = Vl, V^ U V^' = V2 


O3 


Li 




# Wi=A'^3, # W2=Af3 


U w" = Wl, W2 U W2 = W2 
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The lengths of the left subchains are 



Li3 = U^i + Ls ~ L2) , (2.4) 
L12 = ^(Li + L2 — -L3) , 



The structure constant reads 



C^S = y^^^ E Hj^H>H?J(u-|v')(v"*|w')(w"*|u'), (2.5) 

u',v',w' 

where J\fi are the norms of the Bethe states^: 

Ml = (u|u) , M = (v|v) , Ms = (w|w) . (2.6) 

The Hp factors are given by 

H" = S'i(u/,uo) S'i(u2,ui)S'^(ui,Ui)S'^(u2,U2) (2.7) 
Hp = 5i(v'/,vo) S'i(v2,vi) 5^(vi,vi) S'|'(v2, Vg) 

= 5i(w'/, wo) 5i(w^',wi) S>{wi,w';) S>{w2,w'^) . 

with uo = {0^1+^}, vo = {0^2+i| and wo = {0^3+i}_ in the previous formula we used the 
following notations: we denote the scattering factors as 

Sa{Ua,i, Ub,j) = ' _ 2 , a = 1, 2 , (2.8) 

and, given a function F{x,y) and two sets of variables u, v, we define 

F(u,v)= H F{ui,vj), F>(u,v)^ II F{u„vj). (2.9) 



Proportionality factor between ABA and CBA Bethe state: |u)aig = Cu |u)cor is given by 

C„ = [] Y[ ''a,j-Ua,k + i _ ^2.10) 

a=l,2j<it ""'J ""'^^ 

While the formula (2.5) can be explicitly used for a small numbers of magnons, it is not adapted for 
taking the classical limit where the number of magnons is large. The main obstruction for taking the 
classical limit of (2.5) is that the scalar products involved are between off-shell states, and there is no 
closed form expression such as a determinant for this scalar product. In the following sections, we 
restrict our attention to a particular situation where the 3-point function can be written in terms of a 
scalar products of an off-shell state and an on-shell state. 



'in this section all scalar products and norms are understood in the Coordinate Bethe Ansatz normaUzation. 
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3 The su{2) cubic vertex in terms of scalar products 



In preparation for the computation of the type-{2, 3, 3} su{3) 3-point function that we are interested 
in, we review an analogous computation of a type-{2, 2, 2} su{2) 3-point function in [10]. Consider 
the 3-point function of the operators Oi, of lengths Lj, and rapidities Uj with cardinalities Ni, i G 
{1,2,3}. In the following we set ui = u, U2 = v, and U3 = w. In our conventions, Oi consists of 
the fundamental fields {Z, X}, O2 of {Z, X}, and O3 of {Z, X}. 

It is advantageous to generalize the problem slightly by introducing inhomogeneities associated 
with the sites of the three spin chains. Thus the i-th chain is characterized by inhomogeneities 
= {9^\ . . . , 9^1^ }, i = 1, 2, 3. The three sets of inhomogeneities are not independent, because 
the inhomogeneities associated with two sub-chains whose fundamental fields are contracted should 
match. The independent inhomogeneities associated with the contractions between the i-th left sub- 
chain and the j-th right sub-chain are denoted by 0^'^^\ The cardinality of the set 6^'^^^ is Lij. In this 
notation 

0(1) ^ 0(12) u 0(13)^ 0(2) ^ 0(12) u 0(23)^ 0(3) ^ 0(13) ^ 0(23)_ 

The planarity of the (ZZ) contractions between the operators O2 and 03 and the {XX) contrac- 
tions between the operators Oi and O3 selects the component of O3 with N-^ = Lis successive X's 
and L3 — A^3 successive Z's, as in Fig. 1. Consequently, the correlation function is given by the 
product of two factors: 

• The probability to find the component Tr[Z^23^ii3] jn fj^Q gi^^te |w^). 

• The contribution of the remaining contractions can be recast as the scalar product of an on-shell 
state of rapidities u and an off-shell state of rapidities v, in a spin chain of length Li. 

We present below the derivation of the two factors using the language of the six-vertex model. 

3.1 The su(2) Bethe states as six-vertex-model partition functions 

The three type of vertex configurations, a, b, c represented in Fig. 6 have weights 

a(u — z) = -'— , b(u — z) = 1 , c(u — z) = . (3.2) 

u — z u — z 

The rapidities u and z = 6 + i/2 are associated respectively with the horizontal and with the vertical 
lines. These weights are given by the three types of non-zero elements of the L-matrix, which in our 
case coincides with the i?-matrix, 

L{u-z)=l+^—^ , z = e + i/2. (3.3) 
u — z 

Here I is the identity matrix and P is the permutation matrix. 

Consider the expansion of a Bethe vector |u) in the local basis |si, . . . , s^), where Sk G {1, 2}, 

|u) = ^ V.i,...,Si(u)|si,...,SL) • (3.4) 

Sl,...,si,=l,2 

Each of the components V'c,.^ . ^^(u) is a sum over all the possible vertex configurations with on a 
L X N rectangle, with all indices fixed to 1 on the left and the upper boundaries, 2 on the right 
boundary, and free indices equal to si, . . . , s^: on the lower boundary, as shown in Fig. 7. Similarly 
the dual Bethe state is represented by the partition function of the six-vertex model on a rectangle, 
with boundary conditions 1 on the right and the lower boundary, and 2 on the left boundary. 
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2 1— 2 1 1^2 

1 2 

b c 



Figure 6: Graphical representation of the six non-zero el- 
ements of the L -matrix (3.3). The rapidity u is associated 
with the horizontal line, while the rapidity z = 9 + i/2 is 
associated with the vertical line. 



3.2 The su{2) scalar/ inner product in terms of the 6-vertex model 

With the normaUsation of the basis 

( Sl, . . . , SL\ri, ...,rL)= 5si,ri • • • Ss^,rL , (3-5) 

the scalar product of two (in general off-shell) Bethe states 

(v|u)= V's^..^ Ju) Jv)* (3.6) 

si,...,sl=1,2 

is obtained simply by gluing two such partition functions, as shown in Fig. 8, and summing over the 
free indices. In addition to the sesquilinear form (3.6), which is 'the scalar product' , we define a 
bilinear form, which we call 'the inner product' 

(v,u) = (v,u)= Vsi...sl(u) V'st...si(v). (3.7) 

si,...,SL=l,2 

The vertex representation of the inner product (3.7) is obtained by gluing two lattices as the one shown 
in Fig. 7 and summing over the free spin indices. The result is the six-vertex partition function on a 
lattice with indices 1^ on the upper and lower boundaries, and 1^2^ on the left and right boundaries 
in Fig. 8. The symmetry (v, u) = (v, u) of the inner product follows from the symmetry of the 
weights a, 6, c of the vertices in Fig 6 with respect to a rotation by 180 degrees. 
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Figure 7: A six- vertex configurations for the coeffi- 
cient ■0121221211111 (u) °f the Bethe state |u), Eq. (3.4). 
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Figure 8: A six- vertex configurations for the 
inner product ( v, u ) . 
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It follows from the hermitian conjugation properties of the creation (B) and the annihilation (C) 
operators (see the historical note [38]) that for A^-magnon states 

(v|u) = (-l)^ (v*,u), (3.8) 

where the set of rapidities v* is obtained from v by complex conjugation. Since the Hamiltonian 
of the XXX chain is hermitian, the sets of rapidities of the Bethe eigenstates are symmetric under 
complex conjugation. Therefore the normalisation factor Mi = | ( u, u ) | in (2.1) is equal to the 
(squared) norm ( u|u ) of the Bethe eigenstate. 

The structure constant c[^^ is equal, up to the normalization factor, to the cubic vertex made of 
the wave functions of the Bethe states in the representation (3.4), 

C[% = , ^"'"'"^^ = , (3.9) 
\/(u,u) ( V, v) (w, w) 

where the form of the cubic vertex depends on the choice of the three su(2) sectors. In our particular 
case 



(U,V,W) = 2^V'i...iSi^^...si(v) V'.i....^^ 

where the summation indices si, . . . sl^j values 1 and 2. 



2^-2 

^13 



(u) V'2...2 (w) 



(3.10) 
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Figure 9: The su(2) cubic vertex in terms of six-vertex configurations. The shaded area 
factorizes out and can be cut out. Furthermore, the right piece is connected to the rest 
only with type-2 contractions and factorizes out. The lattice splits into two pieces which 
do not talk to each other and can be evaluated separately. 



9 



The cubic vertex ( u, v, w ) can be evaluated using the fact that it gives the partition function of the 
six-vertex model on a lattice obtained by gluing three rectangular lattices with dimensions Li x Ni, 
L2 X and L3 x N3 as shown in Fig. 9. The indices 1 and 2 are identified with Z and X or their 
complex conjugates, depending on the operator under consideration. First we notice that in the part 
of the lattice that has vertical lines labeled by 6^'^'^\ represented by the shaded area in Fig. 9, there 
is only one six-vertex configuration, and therefore its contribution to the cubic vertex factorizes out. 
The factor is a pure phase if the sets v and 0(^3) symmetric under complex conjugation. We will 
assume that this is the case and will ignore this phase factor. Therefore we can delete this part of the 
lattice. 




Figure 10: Schematic representation of the lattice obtained by gluing the rectangular 
lattices corresponding to the states (u|, |v) and |w), with subsequent removal of the 
redundant piece and separating the two non-interacting sub-lattices. After the removal of 
the redundant piece, the states |v) and |w) are no more Bethe eigenstates, because the 
chain are shortened (Li — > L12 and L23)- 

Next, we observe that the sub-lattice associated with the operator O3 factorizes because all lines 
that connect it with the rest of the lattice are of type 2. (This factorisation is obvious in the expression 
(3.10) for the cubic vertex.) These operations are schematically represented in Fig. 10. 

The problem boils down to the calculation of two independent six-vertex partition functions, 
which give the two non-trivial factors in the structure constant. These two factors will be computed 
using the freezing procedure. The freezing procedure for the first factor works as follows. One starts 
from a rectangular lattice corresponding to the scalar product ( v|u ). Both sets of rapidities have car- 
dinality A'^i. The first N2 rapidities v coincide with the rapidities v characterizing the operator ©3, the 
rest A'^s = N1—N2 = L1—L12 of the rapidities v will be denoted by wtvj+i = z^-^^^+ii ■ ■ ■ i^Ni = -Sli. 
or symbolically, v = v U z. 

3.3 The su{2) freezing procedure 

If we adjust the rapidity of the last magnon to the value of the last inhomogeneity, z^-^ = Ol-^ — 
then the vertex at the low right corner is necessarily of type c. Then the only possibility for the rest of 
the vertices on the last row and the last column is that they are type b. This is what we call "freezing". 
Hence last row and the last column form a hooked index line carrying the index 2, as shown in Fig. 1 1 , 
left. This procedure is repeated N->, times, the rapidities of the lowest N'i rows fixed to z = 0(13) _^/2. 
The result is that the rightmost N'j, indices below the lowest horizontal n-line are fixed to the value 2. 
After removing the frozen part of the lattice, shaded in blue in the figure, we obtain that the first factor 
in the cubic vertex equals the scalar product ( v U z, u ). The contribution of the frozen part of the 
lattice is the product of all c-vertices on the diagonal, which equals (— l)^i^, a factor we will ignore. 

In a similar way we compute the second factor . The freezing procedure is shown in Fig. 1 1, right. 
We start with a scalar product ( vi^, z ) for a chain of length = + L23. We freeze the rapidities of 
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Figure 1 1 : The freezing procedure for tlie two factors in C123 



the bra state to z = 0^^^^ — i/2. The frozen area (shaded in blue) gives a contribution, which is a pure 
phase if both sets w and 6^"^^^ are symmetric under complex conjugation. We will assume that this is 
the case and will ignore this factor. The rest of the lattice gives the second factor in the expression for 
the cubic vertex. We find 

(u,v,w) = (vUz,u)0(i, (w,z)0(3) (3.11) 

up to a factor which takes into account the contribution of the deleted and added pieces of the lattice. 
This factor is a pure phase, since the set w is symmetric under complex conjugation, and can be 
ignored. 

4 The su{3) cubic vertex in terms of scalar products 

As before, we consider that the three operators, Oi, O2, O3, are described by three sets of rapidities 
u = {ui, U2}, V = {vi, V2} and w = {wi, W2} with cardinalities respectively A^i + Mi, A^2 + M2 
and + M3. In the configuration we are considering, W2 = 0, since O3 is an su{2) operator. We 
refer to Section 6 for the equations obeyed by these rapidities. 

Again, we have two types of contributions to the correlation function: 

• the contribution of the (ZZ) contractions between the operators O2 and O3, through the factor 
(zi,wi), withzi = 0(13) _^/2and #zi = N^, 

• the remaining contractions, which can be recast as the inner product (v U z, u) between an 
on-shell vector of a spin chain with length Li and rapidities u = {ui, U2} and an off-shell state 
with the same length and rapidities (v = {vi, V2}) U (z = {zi,Z2}), with zi = 9^^^^ — i/2 
and Z2 = 0(13) - i. 
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Below we evaluate, using the freezing argument, the {2, 3, 3} type structure constant (Fig. 1), 



r 



(0) 
123 



(u, V,Wi 



si[{3) 



su(3) su(3) , 

(u,u) (v,v) (wi,wi 



(4.1) 



We will show that the corresponding cubic vertex is given by 

(U,V,W) = (VUZ,U)„{1) (wi,Zi)0(3) 



(4.2) 



with 



(4.3) 



Here ( , ) denotes, as before, the su{2) inner product, and ( 



denotes the sii(3) inner product. 



4.1 The su{Z) Bethe states in terms of the 15-vertex-model 

In order to generalize the freezing procedure to su{2>), let us first show how to represent the compo- 
nents of the su(3) Bethe vectors in terms of configurations of a 15-vertex model shown in Fig. 12. 
The vertices are similar to those from Fig. 6, with the difference that the indices carried by the lines 
can be now 1 , 2 or 3. We represent them graphically by thin, red and black lines, respectively. The 
weights are identical to those from equation (3.2), depending on whether the indices carried by the 
lines are equal or different. 



2 I 2 

2 



3 I 3 



1 

a 




2 1 

2— i 1 1 ^2 

1 2 

3 1 

3— L-1 1— p 3 



2 3 

2 I 3 3 I 2 

3 2 



Figure 12: Graphical representation of the 15 non-zero elements of the su(3) 
L -matrix, Eq. (3.3). The rapidities u and z ~ 6 + i/2 are associated with the 
horizontal and the vertical lines, respectively. 



The Bethe vector | u) is given by the expansion 



u 



si,...,sl=1 



^si,...,Si(u) |S1, 



(4.4) 



where V'si,...,sl (") ^ ^^^^ "^^e possible 15-vertex configurations on a rectangular lattice with 
Li + A^i vertical lines and A'^i + Mi horizontal lines, with the free spin indices equal to si, . . . , s^. 
An example for such a vertex configuration is given in Fig. 13. The first Li vertical lines carry 
rapidities 9^^^ + i/2,..., 9^^^ + i/2 and spin indices 1 on the top, which correspond to the vacuum 
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|0) = = |ll...l). The right vertical lines carry rapidities till +i ... uiAT^+z and have index 
2 on the top. At the bottom, the first Li indices are free, and the last A^i ones are fixed to 1. The lower 
A^i horizontal line correspond to the first-level magnons and carry rapidities ui^i, . . . ,ui^Ni- The 
higher Mi horizontal lines represent the second-level magnons with rapidities U2,i, . . . , U2,Mi- Due 
to the particular spin and rapidity choices, the shaded regions are frozen to the particular configuration 
shown in the Figure. This diagram is equivalent to (a special case of^) the one used by Reshetikhin 
in [39]. 



n,+ i/2 \ 2 ■ 



Ui 



^<"+«/2 



111111111111 



2 2 2 2 



131231311111 



1111 



Figure 13: A configuration contributing to the coefficient V'i3i23i3iiiii(u) of a su(3) 
Bethe vector (4.4) with L = 12, TV = 4 and M = 3 . 



The structure constant factorizes as in the su{2) case. The two factors can be cast in the form of 
scalar products of an on-shell and an off-shell Bethe states by applying the freezing procedure. 



4.2 The su{3) freezing procedure 

Consider the scalar product of two su{3) states of the first chain, (v, u)=({vi, V2}, {ui, U2}), as 
represented in figure 14. Our purpose is to freeze the rightmost N3 indices to the value 3, as imposed 
by the planarity of the contractions in the three point function. This can be done by setting the last N3 
rapidities of first level magnons at their freezing values, 

{ii,7V2+i = ^i'^^ - . . . , vi,N^ = 9^N^^ - i/2. (4.5) 

This will insure that the corresponding frozen region contains only red and black lines propagating 
from the top to the bottom of the diagram. The number of black or red lines is not fixed by the 
freezing, only their sum is fixed. In order to force all the lines in the frozen region to be black, we 
have to apply once again the freezing procedure to the magnons of second level, by fixing 

V2,M2+1 = ^i'^^ - • • • ^2,Mi = 6*^3^^ - i- (4.6) 
The remaining magnons are set to the corresponding values in the state | v), 

■I'l.l = -Vl.l, ■ • • , V2,M2 = V2,M2- (4-7) 

This gives us the first factor of the expression (4.2) for the cubic vertex. The second factor, ( zi, wi ), 
is the same as in the su{2) case. 

'^In [39], there are had two momentum-carrying nodes, while our spin chain has only one momentum-carrying node. 
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ui 
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1 



• Vi 



}e(i3)-i/2 



2 } V2+1/2 
2}0(i3)-i/2 



2222 111111111111 

vi+i 6)(i3)+i/2 6>(i2)+i/2 e(")_^i/2 

Figure 14: The inner product (v, u) and the freezing to (v U 0^^^^ — i/2, u) 



5 The su{2) structure constant in terms of the -functional 

In the su{2) case, the rapidities of the on-shell Bethe states satisfy the Bethe equations 

fl^^l^p^ fl^^l^lllll ^ j = l,...,N. (5.1) 

The Bethe equations (5.1) follow from the requirement that the eigenvalue T{u) of the T-matrix, 

T{u) = ^"'(^^ + ^"""(^^ ^ e^f i(-) + e'Piin) (5.2) 

Qe (u) Qu{u) Qu{u) 

has vanishing residues at the Bethe roots. ^ Here Qu and Qg are the Baxter polynomials 

N L 

Qu{u) = ll{u-Uj), Q0{u) = '[l{u-ei), (5.3) 

j=i 1=1 

and 

Q^iu) = Q{u±i/2), Q±±(m) = Q(n±i), gM(n) = Q(u + m/2). (5.4) 



'with the normalisation (3.3) of the L-matrix, T{u) is not a polynomial, but has poles at u 
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The Bethe equations read 

g2ipuK)^_l^ j = l,...,N, (5.5) 
where the pseudomomentum known also as the counting function, is defined modulo tt by 

Qq Qu 

The T-matrix (5.2) is normaUzed according to the vertex representation with weights (3.2). The 
eigenvalues of the diagonal elements of the monodromy matrix on the vacuum, A{u) and D{u), are 
given by 

L 

D{u) = e^fe^") = Ylb{u-9i- i/2) = 1, 
A{u) = e^^'eW = a(ix - - i/2) = ^^^""^ 



(^) 



The inner product of an on-shell Bethe state u and an off-shell state v is given in a determinant 
form by [40] [41] 

N 

(v,u) = []Ae(t;,)^u,v, (5.8) 



where 



is the Slavnov determinant. The Slavnov kernel Q,{u, v) is defined by^ 



n{u, v) = t{u -v)- e^^P"^") t{v - u) , t{u) = ^ - . (5.10) 

The Slavnov determinant cannot be directly evaluated in the classical hmit. This can be done 
using the representation in terms of the =(2/-functional introduced in [29]. The Slavnov determinant 
-^u,v was expressed in terms of the -functional first for the Umit u — > oo [29] and then in the 
general case [32] [31]. Later a more compact expression was found in [37]. 

The j2/-functional, whose properties are listed in Appendix C, is defined for any function f(u) 
and any set of points in the complex plane u = {ui , . . . , uat} as follows, 

det,fe(uJ-^-/(n,)(u,±zM 

<[f] = — rr^^ ■ ^^-^^^ 



^In order to simplify tlie formulas, liere (as well as in [32] [31] [37]) we use a different normalisation for the Slavnov 
matrix than in [41]. In our conventions, the Slavnov kernel depends only on the pseudomomentum pu = Pu — Pu- 



15 



The Slavnov determinant is expressed in terms of this functional as [37] ' 

■^u,v = (-l)^<+uv[%]- (5.12) 

Let us express the scalar products in the expression for the cubic vertex in terms of the s^- 
functional. We find for the two inner products in (4.2) 

(vuz,u)^a) = (-i)A^^nA(.,)n^(%Xuu[S^], (5.13) 

j=l j=l *^6>(i2) 

Ns Q+ 

(w,z)^(3) = (-l)^3j]A(z,-Xr[^]. (5.14) 

j=l ^0(23) 

Proof: Using the properties of the ^ -functional (Appendix C), we transform 

"^uzuuIt^J - "^uzuuiT^^Ti^J --^uuIt^^J^ (5.15) 

Wgil) ^0(13) ^0(12) ^0(12) 

^+ = ^+ r Qe(13) Qe(23) , +. Qg(23) , 

Ignoring the factors that are pure phases, we find for the cubic vertex 

(u,v,w) =<uu[7;^]<r[7:|q— ]■ (5.17) 

^0(12) ^0(23) 

The second factor has been evaluated in a different ways in [10] and [31], where it was used that it 
equal to a partial domain wall partition function of the six- vertex model. 

The norm ( u, u ) is most easily computed by taking the expression for the inner product ( v, u ), 
Eq. (5.12), in the limit v ^ u. 

6 The su{3) structure constant in terms of the -functional 

A generic Bethe state |u) in an su{3) sector is characterized by the rapidities u = {ui, U2} and the 
inhomogeneity parameters 6 associated with the momentum-carrying node (1), where 

Ui = {uij,...,Ui^n}, U2 = {u2,1, . . . ,U2,m} , 9 = {9i, . . . , Ol} ■ (6.1) 

The rapidities satisfy the nested Bethe wave functions for the 5^(3) R-matrix given by (3.3): 

L zi , 1 • N , ■ M 1 

UlJ -dl + 1^1 _ -p-r UlJ - Ul^n + I -TT - ^2,1 



2^ 



TT "Ij —^1^2'' _ TT "Ij ~ "l,™ ' TT 

f=l Ul,3 „=i - "1," - i iii - U2,m + \i 

M . . N 1 . 

'The expression (5.12) for the Slavnov determinant depends on the ensemble of the rapidities u and v in a completely 
synmietric way. This remarkable symmetry of this expression follows from the fact that, due to the global su{2) symmetry, 
the annihilation operators with rapidities v can be replaced by creation operators with the same rapidities and a global 
raising operator [37]. See also the exercises of the 3rd day of the 4th Mathematica School (http://msstp.org/?q=node/272). 
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The Bethe equations (6.2) follow from the requirement that that the su{3) T-matrix in the funda- 
mental representation, 



T{u) = e^J'"(") + e'*'"(") + e*^'"^") , (6.3) 

has vanishing residues at the Bethe roots. For given distribution of the roots ui and U2, the pseudo- 
momenta are determined modulo 27r by^ 

Qui Qq Qui Qu2 Qu2 

see e.g. [42]. In terms of the three pseudomomenta, the Bethe equations (5.5) read 

e^PuW-ipSW = _1 if^Gui; 
^Pl{z)-ipl{z) if ^ - i/2 G U2 . 



(6.5) 



It is convenient to introduce the functions Pl^{z) and Pl^{z), associated with the two nodes of the 
Dynkin graph of su(3), and related to the quasimomenta i = 1, 2, 3, by 

^u(^) = Pu(^) - Pu(^), ^u(^) = pI{z + V2) - Pi{z + i/2). (6.6) 

In terms of these functions, which we will also call pseudomomenta, the Bethe equations take the 
more standard form 



e*-PuW = _l^ ifzGUa (a = 1,2). (6.7) 
The functions Pi and P2 can be expressed in terms of the su(3) Cartan matrix {M^b} = (_i "2) 

6=1,2 Quj, 



jpa 

e " 



Q^J U^^r « = 1'2. (6.8) 

^0 ^ 6=1.2 



Let us stress that the values of the local conserved charges are determined only by the level- 1 roots 
ui. The duaUty transformations change the level-2 roots U2, but leave invariant the level- 1 roots ui, 
which carry the physical information [42]. 

The norm of an on-shell Bethe state The squared norm of an on-shell Bethe state has been computed 
for the case of su(3) by Reshetikhin^ [39] and is expressed as the determinant of the matrix of the 
derivatives of the two quasi-momenta: 



(u,u) = Cu det du^ .P^{ub,k] 



(6.9) 



where the determinant is with respect to the double indices A = {a,j} and B = {6, k}. The nor- 
maUzationn factor Cu is given by (2.2). The matrix of the derivatives of the two quasimomenta is 
explicitly 



dP"'(z] 



dz 



, (6.10) 



Here we used the conventions of Eqs. (5.3) and (5.4). 
'a conjecture for su{n) is proposed by EGSV in [9]. 

'"Here it is assumed that the set of the Bethe roots is symmetric under complex conjugation. 
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where 

tab{u) = - . (6.11) 

Instead of taking the derivatives, we will compute the norm as the limit of the determinant depending 
on two sets of rapidities, u and v, which has the hmit (6.10) when v — > u. We define the (iV + M) x 
(A^ + M) square matrix ^a},{uj,Vk), with 

9.ah{u,v) = tabiu-v)-e'^-^''Uab{-U + v). (6.12) 

The expression for the norm, which we are going to evaluate in the classical limit, is 

( u, u ) = Cu lim det [0.ab{ua,j , Vb,k)] ■ (6.13) 

6.1 The inner product (u, v) in the limit U2 — >^ oo 

Unhke the su{2) case, the inner product of an on-shell Bethe state with an of-shell Bethe state is not 
generically a determinant. Determinant representations exist in some particular cases [36,43,44]. We 
will use the determinant expression obtained by Wheeler [36], when the rapidities of the second type 
of magnons of the Bethe eigenstate are sent to infinity. We assume that M is odd; then one can send 
to infinity the U2 roots one by one. As a result the second level Bethe equations become trivial and the 
first level Bethe equations take the same form as for su{2). The inner product ( u, v )g factorizes 
into two su{2) inner products [36] 



^lim (V, u)-^^' = det - - + 

U2^oo jk \ Qvi{,V2,j)J 



'J \ Qe{n,k) Qui {vi,k) 



= (ui,vi)r(') (oo,V2)tf). (6.14) 
Using (C.13), we write (6.14) in the form 

lim (v,u)e = <uvj^] <[^] • (6.15) 

U2->oo (y^^ 

The su{3) structure constant in terms of the j2/-functional. Combining (5.13), (5.14), (6.14) in 
Eq. (4.2), we get 

(u,v,w)=<[§f^]<,,J§^] (6.16) 
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7 The semi-classical limit of the su(3) 3-point function 



7.1 The Sutherland Umit 

The classical , or thermodynamical limit is attained for long spin chains (L ^ 1) with macroscopically 
many excitations N,M ^ L, and in the low energy regime {5E ~ [1 1, 12,45]. Such spin chains 
correspond to "heavy" operators, which are traces of products of many SYM fields. In this Umit 
the roots scale as m„j ~ L. In the condensed matter literature the classical limit has been studied 
by Sutherland [46] and by Dhar and Shastry [47], and is known as Sutherland scahng limit. In the 
classical limit the roots are organized in several macroscopic strings, which condense into cuts in the 
complex rapidity plane. The three quasimomenta pi,P2,P3 become the three branches of the same 
meromorphic function. The three sheets of the corresponding Riemann surface are joined among 
themselves along the cuts defined by the long Bethe strings. In the classical hmit, the Bethe state is 
characterised by the resolvents 

GuAz) = dz\ogQu{z), Gu^iz) = dzlogQu^iz), (7.1) 
as well as the resolvent for the inhomogeneities 

Gg{u) = dulogQeiu). (7.2) 

The two resolvents, Gu^ and Gu2, can be expressed in terms of the three quasimomenta Pu,Pu and 
Py, which become the three branches of a single meromorphic function on the tri-foliated Riemann 
surface, 

Pu = - Ge (mod 27r), 

p2=Gu2-Gui (mod27r), (7.3) 
P'L = -Gu2 (mod27r). 



or 



Pi = 2Gu, - - Gg (mod 27r), 
P^ = 2Gu2-Gui (mod27r). 



Let Cfj be the cuts joining the i-th and the j-th sheets. Then the Bethe equations (6.5) become 
boundary conditions on these cuts, depending on the mode numbers nfy. 

where ]/) denotes the half-sum of the values of the function p on both sides of the cut. 



7.2 Stacks 

In addition, there is the possibility of configurations called stacks (bound states of rapidities associated 
with different nodes [48]), which represent pairs of roots belonging to the nodes 1 and 2 and at distance 
0(1) from each other [42,45,49]. We can have macroscopic strings of stacks, which in the classical 
limit become two cuts that merge into one cut. Since the roots that form the string of stacks belong to 
two different nodes, they correspond to a cut type 1-2 and a cut type 2-3, where we understand that the 
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cut of type i-j joins the i-th and the j-th sheets of the Riemann surface. The resuh of merging of the 
two cuts is a cut of the type 1-3. Therefore, in order to have a description of the generic Bethe state in 
the classical limit, we must assume also the existence of cuts of type 1-3. The boundary condition on 
these cuts is obtained by taking the limit of (7.5) and has the form 

27r<3=;^i-^, ^eCfg. (7.6) 

The bosonic duality transformations [42] in the classical limit corresponds simply to the exchange of 
the Riemann sheets 2 and 3. 

7.3 The semi-classical norm 

The determinant (6.13) can be computed in the classical limit under the assumption that there are only 
1-2 and 2-3 type cuts, which are separated at macroscopic distance ~ L. With this assumption, the 
off diagonal elements of Q{u, ~ (u — v)~^ ~ and the only matrix elements of order one are 
those in a strip of width ~ 1/ VX along the diagonal. As a consequence, the non-diagonal blocks do 
not contribute in the classical limit and the determinant is simply the product of the determinants of 
the diagonal blocks, 

(u,u),„(3) ~ det det[n22{u2,j,U2,k)] 

Let us evaluate the norm assuming that there there are no cuts relating the first and the third sheet of 
the Riemann surface. We can use the expression for the classical Umit of the norm in the su{2) sector 
(Appendix (C)): 

(u,v), = log<+,[^]= f gLi2[/(z)e^«"(^)+^«v(.)-«G«(.)]+^(l^g^)^ 

with 

Gu(z) = d,logQ^{z), Gg{z) = d,\ogQB{z). (7.9) 
The norm of the classical Bethe state is then^^ 

log(u|u)=/ ^Li2(e2^Gui(^MC?u2(z)-iGe(z)-) ^ / ^Li2(e2^Gu2(^MGui(z)V g ^^^ 

/Cui 27r Jc^^ 27r 



*'We conjecture that in the most general case, when some of the roots can form bound states ("stacks"), this logarithm of 
the norm is given by 

log(u|u> = ^f |iLi2(e-SW-^^gW), 
where Cij {i, j = 1,2,3) denote the contour (or contours) surrounding the cuts between the i-th and the j-th sheets. 
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7.4 Semi-classical limit of the structure constant 



Substituting (7.8) in (6. 16) we find for the structure constant in the classical Umit 



+ / ^ Li2 (e^G.,(z)-iG., (.)) ^ / ^ Li2 (e'^-i (^)-^Gei^s) (-)) 



(7.12) 



- |log(u|u) - 5log(v|v) - ^log(w|w) . 



The last Une is evaluated in the classical limit according to (7.1 1). 

8 Conclusions and outlook 

We have analyzed the tree-level 3-point functions of single-trace operators of the planar Af = 4 SYM 
theory in the su{S) sector. Each of the three operators is an eigenstate of the dilatation operator, and 
it is characterized by a set of charges (angular momenta) and a set of rapidities. We have classified 
the possible configurations of Wick contractions and given the general expression of the 3-point func- 
tions in terms of rapidities associated to each operator. This expression, obtained using the tailoring 
technique of EGSV [9], is not adapted for taking the classical Umit. In some particular situation, 
when one of the operators belongs to an su{2) sector, we are able to express the 3-point function 
using the alternative method of freezing proposed in [10]. By further specializing the second group 
of rapidities corresponding to one of the operators, we can use a result of [36] to express the scalar 
products as determinants. Finally, the semi-classical limit of the determinants can be taken using the 
results from [3 1,32]. The simplest classical operator from the su{3) sector are the three-spin solution 
obtained by C. Kristjansen [50]. 

There are two obvious directions to explore. First, one can try to evaluate the general sum over 
partitions in (2.5) quasi-classically. One can either try to perform a quasiclssical evaluation of the 
Korepin sum over partitions for the scalar product of two off-shell Bethe states, in the spirit of [29], or 
refine the coherent state approximation [25, 51, 52]. Another direction is to explore the non-compact 
sector of the theory. There are several recent papers which are relevant for that, [53-59, 59-62]. 
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A The nested coordinate Bethe Ansatz 

As an introduction to the 'tailoring' procedure in the context of su{3) 3-point function, this appendix 
is a brief introduction to the nested coordinate Bethe anstaz of the su(3) spin chain. Let us consider 
an su{K) spin chain of length L. The treatment follows closely [63]. The hamiltonian reads 




(A.1) 



n=l 
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At each site of the spin chain, there is a spin with K different polarizations. The Hilbert space of 
the spin chain is "H = (C^)®^. In the Hamiltonian (A.l) In,n+i is the identity operator in the space 
C^_,_]^ and Pn,n+i is the permutation operator in ® '^n+i- Dynkin diagrams of the Lie 
algebra give a convenient way to label the excitations. For the su{K) algebra, the Dynkin diagram is 
given by Fig.(15). 

cwx^ -o 

12 3 K-l 

Figure 15: The Dynkin diagram of su{K) algebra. 



Following Bethe, one looks for the eigenstates of the Hamiltonian (A.1) in the form 

I*) = X] V'ni,...,nK-i(uir-- ,Ui^_i)|ni,--- ,nx_i) , (A.2) 

positions 

where Uq is the set of rapidities of node a and ria labels the positions of the excitation of node a 
{a = 1, - ■ ■ ,K — \). The summation is taken over all possible positions of excitations. The nested 
ket state |ni, • • • , riK-i) is constructed in the following steps: 

1 . Start with an initial value of length L; 

2. Create A^i excitations of node 1 at positions ni. The A^i excitations form a reduced inhomoge- 
neous spin chain of length Ni, 

3. Create N2 excitations of node 2 at positions 112 in the reduced spin chain. One should have 
N2 < Ni. The excitations of node 2 again form a reduced inhomogeneous spin chain of length 

4. Repeat the above steps for all the K — 1 nodes. 

The procedure for the su{3) case is explained in Fig.(16). We would like to stress the important 
feature that the excitations of node a should only be created from the excitations of node a — 1 (the 
reduced spin chain). Therefore we have Na < Na-i. The positions should obey 

1 < m,i < ■ ■ • < ni^Ni < L, 

1 < na,i < ■ ■ ■ < na,Na < Na-i, foT 2 < a < r . 

The beauty of this method is that at each step, the operation is simple and the same: to create one kind 
of excitation in a reduced spin chain of length Na. Here we see the nested feature of this method. 

The wave function ip{ui, • • • , uk-i) is constructed by nested Bethe ansatz. It is given in terms 
of a series of wave functions ijja, a = 2, • • • ,K — l.The Bethe wave function reads 

V'lui,--- ,Ui,_i) = ^Ai(Pi)n MPi) , (A.3) 

where the wave functions a = 2, • • • , -fC — 1 are given by 



22 



Initial vacuum of length L 



Create Ni excitations at positions {nij} 

-O^-O^OO-O^OO-O^-O^O-O^O- 

Create A'2 excitations at positions {n2j} in the reduced chain of length A^i 

o o o^^^ — o — 

Then the final state {"2,^}) is represented by 



Figure 16: Nesting procedure for su(3). In our example, K = 3, L = 12, {nij} = {1,2,5,8,9,11} and 
{»^2j} = {1, 3, 4, 5, 6}. The final result reads |2, 1, 0, 0, 2, 0, 0, 1, 2, 0, 2, 0), where denotes the vacuum state, 
1 denotes the first excited state and 2 denotes the second excited state. In our ket notation, it is written as 

|{1,2,5,8,9,11},{1,3,4,5,6}). 



with Pa the permutation {Pa,i, Pa,Na) of (1, • • • , Na). We define V'k(Pk-i) = 1- Ma,b is the 
Cartan matrix of the Lie algebra, for su{K) being 

Ma,b = ^Sa,b - Sa-l,b - Sa+l,b ■ (A.5) 

Our choice of normalization is 

Aa{l,2--- ,Na) = l (A.6) 

and the coefficients Aa obey the relation 



Here we used the following definition, 



S2iUa,i,Ua,j) ■ (A.7) 



Sa{Ua,i, Ub,j) = '_l , £7 = 1,2. (A.8) 

In order that (A.2) is an eigenstate of the spin chain Hamiltonian, the rapidities should satisfy the 
Bethe ansatz equations : 

/ UaJ + Va£\ ^ ^ TT TT ^a,j " Ub,k + 



where Va are the Dynkin labels. We consider the fundamental representation in this paper, where 

Va = Sa,l. 

B The su{3) tailoring prescription 

We consider the operators in sn(3) sector with definite one-loop anomalous dimensions. In the spin 
chain language, these operators are represented by the Bethe eigenstates of the su{3) spin chain. 
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One first write the Bethe state as the entangled state of two subchain states. This operation is called 
"cutting". After cutting operation, each subchain state also takes the form of Bethe states. In order 
to perform Wick contraction, one needs to "flip" one of the subchains. Flipping is an operation that 
takes a ket state into the corresponding bra state with the same wave function. One can flip either the 
left or the right subchain. In this paper, we always flip the right subchain. The last step is to calculate 
the scalar product of Bethe states, this is called the "gluing" operation. 

Consider a generic su{2) Bethe state |u) of a spin chain with length L. We define the first I sites 
from the left to be the left subchain and the rest L — I sites to be the right subchain. |u) can be written 
as an entangled state of the subchains 



min{iV,«} 

|„)= V 



E 



^ ^p{u)\nl,■■■ ,nk)i^\nk+i-l,--- ,riN-l)r (B-1) 



k=0 l<ni< --<nfc</ l<nk+i<---<L 

where k is the number of magnons in the left subchain. Note that one needs to re-label the positions 
of the magnons in the right subchain, see Fig. (17). Since we have two subchains, the magnons can 



O O O O' 

12 3 4 
left subchain /""^ 



■o o o o 

5 6 7 8 
right subchain 



o o o o o o o 

1 2 3 4 1 2 3 4 



Figure 17: The cutting process for a spin chain with L = 8. We take I = 4 and L — I = 4. The sites in the 
right subchain are originally labeled by 5, 6, 7, 8 while after cutting they are labeled by 1, 2, 3, 4. 



either be in the left subchain or in the right subchain. After cutting a Bethe eigenstate, the resulting 
two subchains states still take the form of Bethe state. Hence the two subchain states have their own 
Bethe wave functions ipi{u') and il)r{vi!'), where u' and u" is a partition of u, 

u' U u" = u, u' n u" = . (B.2) 

In general, V'(u) = i?(u', u") ^/'/(u') i/^rW) where H{u', u") is a partition-dependent factor which 
shall be called i?-factor from now on. Formally, the cutting of a Bethe state can be written as 

|u) = H{u', u") |u') ® |u") . (B.3) 

a 

Usually the expression for the i^-factor is long. In order to make the expression more compact, we 
introduce some short-hand notations. Given a function F{x,y) and two sets of variables u, v, we 
define 

F(u,v)= H F{ui,Vj), F>(u,v)= J] F{ui,vj). (B.4) 

■UjGu, Vj €v 

For a constant c, we define 

F(u, c) = n F{ui, c), F{c, v) = n ^(c, Vi) ■ (B.5) 
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With the notations from (A.8) the iJ-factor for the su{2) spin chain is given by 

H{u', u") = Si{u", Oy S^{u', u") . (B.6) 

The cutting operation can be generaUzed to su{K) Bethe state. Let us denote the nested Bethe state 
by |u) where u = (ui, • • • , uk-i). We have 

|u) = H{u', u") |u') |u") (B.7) 

u' 

with the i^-factor 

K-l 

H{n',u") = n S>{n'^,u';) (B.8) 

n=l 

where Uq is defined as Ug = {O'} with I the length of left subchain. In order to perform Wick 
contraction, we need to "flip" the right subchain from a ket state into a bra state. The flipping operation 
is different from Hermitian conjugate. Given a state 

lip) = e'^\XZXZZ) (B.9) 

the Hermitian conjugate and flipping (denoted by superscript JO lead to 

(|^))t = e-'^{XZXZZ\ 
(1^,))^ = e+'^{ZZXZX\ . 

For an 5-^(2) Bethe state |u), the flipped state is proportional to the hermitian conjugate of the Bethe 
state |u*) 

(|u))^ = (u*|F(u) (B.IO) 

where u* is the complex conjugate of u and we call the proportionality F{u) the F-factor. For an 
su{2) Bethe state, the F-factor reads 

F(u) = 5i(u, 0)^+1 52>(u,u) (B.ll) 

where L is the length of the spin chain. For an su{K) Bethe state |i3), the F-factor is given by 

K-l 

F{U) = n S,{Un,Un-l)S>{Un,Un) (B.12) 

71=1 

where uq is defined by uq = {0^+^}. From now on, by "tailor" (denoted by T) a (nested) Bethe state, 
we mean first cut the state and then flip the right subchain state. We define the product of the corre- 
sponding i^-factor and the F-factor to be the H^^ factor. For an su{K) Bethe state |ui , • • • , u^-i) 

(|u)f = ^H^|u')®(u"*| (B.13) 



where 



and the i^jr-factor reads 



u' U u" = u, u' n u" = (B.14) 



K-l 



= n Sl{ulUn-l)S>{Un,K) . (B.15) 



n=l 
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C The functionals 



C.l Definition 

For any set of points u = {uj]^^^ in the complex plane and for any complex function ^{z), we 
define the pair of functionals [/], which are completely symmetric polynomials of degree N of 
the variables f{ui), ... , /(ujv). The functional is defined as a sum of monomials labeled by all 
possible partitions of the set u into two disjoint subsets u' and u", with u' U u" = u, 

= T. n (-/("Oi n <c.i, 

u'uu"=u u'eu' u"eu" 

The expansion (C.l) was thoroughly studied by Gromov, Sever and Vieira [29]. The expansion (C.l) 
is summed up by the following operator expression, 

<[/]=<^ [/]•!, (C.2) 
where the different operator si/^ [u] is defined as 

" j=l j<k 

The operator functional [u] is formally obtained from the c-functional [u] as 

^u[f] = <[fe'^]- (C.4) 

C.2 Properties 

It was found in [29] that for constant function f{u) = k, the expansion (C.l) does not depend on the 
positions of the rapidities u and the functional [/] is given in this case by 

<±M = (1-Ac)^. (C.5) 

The functional [/] have other remarkable properties [31, 32]: 

1) Determinant formula, Eq. (5.11), which is obtained by substituting Au = detjk{u^~^) in the 
definition (C.3). 

2) Functional relations between s^'^ and s^~: 

N 

i=i 

<[f] = <[-^f]- (C.7) 



3) Reduction formula: 



<uz[/7^] = <[/]• (C.8) 
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4) Factorization property: 

<bv [/] = <^ /] • <^[^ f] ■ (C.9) 
Aplying A'^ times the factorisation property, one obtains the representation 



Relation to the Slavnov determinant. The determinant of the kernel n{u, v) defined as 

nf{u,v) = t{u -v)- ^{u) fiv) ti-u + v), (C.ll) 
where the set u satisfies the "on-shell" condition 



Qu (Uj) 



f{uj) = -1 for Uj e u, (C.12) 



is evaluated as [37] 



C.3 Classical limit. 



In the classical limit, the Bethe roots condensate in one or several disjoint cuts. Let Ck be a contour 
encircling the fc-th cut anticloskwise and leaving outside all other singularities of / and Gu- The 
filling fraction of the A;-th cut is 

ak = ^jGMdu. (C.14) 

We consider the limit L — oo with all finite. Then the leading, linear in L, term of log is 
given by the contour integral 

log<±[/] ±y"^Li2(/(u)e±^«"W), Cu = ULA. (C.15) 

Cu 

While there is not yet a rigorous proof of this formula, it has passed a number of analytical and 
numerical checks. A heuristic derivation of (C.15) for f{u) = ^ ^"^y^ ) > was presented in [29]. 
When f{u) = k, it was shown in [29] that the quasiclassical formula (C.15) gives the exact answer 
(C.5). Moreover, the integral (C.15) satisfies the functional equation (C.6) thanks to the functional 
equation for the dilogarithm, 

Li2(-) = -U2{co) -^-\ \og\-uj). (C.16) 
UJ o z 
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